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Abstract 

In this note, the Sherrington Kirkpatrick model of interacting spins is under consideration. 
In the high temperature region, we give an asymptotic expansion for the expected value of 
some genereral polynomial of the overlap of the system when the size N grows to infinity. 
Some of the coefficients obtained are shown to be vanishing, while the procedure to get the 
nontrivial ones has to be performed by a computer program, due to the great amount of 
computation involved. 
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1 Introduction 

This note is concerned with the usual Sherrington-Kirkpatrick model of interacting spins, that can 
be briefly described as follows: for N >1, set Ev = {—1,1}'^, called the space of configurations. 
For a given a = (ui,..., a^) G E^v, the quantity Ci will stand for the value of the spin of the 
configuration. We will suppose that the spins randomly interact according to a mean field model 
represented by a Hamiltonian given by 

where {gij;l < i < j < N} is a family of independent standard Gaussian random variables, 
defined on a probability space (0,.F, P). For a given positive parameter f3 (that stands for the 
inverse of the temperature of the system), set then 
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We will consider here the random Gibbs measure Gn on E^v, that gives the highest weight to those 
configurations minimizing the energy i?Ar, and whose density with respect to the uniform measure 
is given by exp(—/3i/7v(o’)) for all a G Ejv- For any fV, n > 1 and any function / : —» R, 

we will then set (/) for the average of / with respect to dG®", i.e. 


(/) = ^w" X! /(o-\...,cr")exp 




Kr. 


and the quantity E(/) will be denoted by 

Parisi’s analysis of the SK model (see Q for a complete discussion) showed that a family of 
functions of particular interest, summarizing a great amount of information about the whole spin 
system, is the family {Rij'; 1 < I < I' < oo}, where 

i<N 


and , cr* are understood as two independent configurations under the measure dG®^. The 
quantity Ri^i> (which will be the main object of study of the present paper) is called overlap 
between the configurations cr* and , and can be easily related to the number of spins of cr* and 
which are in the same direction. 

While the efforts of Physicits generally aim at the understanding of the low temperature regime 
of the SK model, and of the related phenomenon of non self averaging (see e.g. D and !§), the high 
temperature domain (i.e. /3 G [0,1)) is still of mathematical interest: indeed, in that case, a rather 
complete description of the limiting behavior of the system when N oo can be given, generally 
amounting to some asymptotic theorems concerning and Ri^i'. Among those results, we can 
cite concerning the L^-limit of ^ log(ZAr), as well as §, i and 1^, that provide some Central 
Limit Theorems for Zn and Ri, 2 - However, the possibilities of getting some sharp estimates for 
the above quantities have been greatly enhanced by the rigorous introduction and proof of the 
cavity method, summarized in , leading among other results to a number of interesting CLT for 
the family {Ri^i'; 1 < I < I' < cxd }. 

In this paper, we would like to go one step further into that line of investigation, and try to give 
a complete expansion (in terms of power series of -^) for the quantity 2 )) more generally, 
for any quantity of the form 

^ ^ ( 1 ) 

for an arbitrary sequence of indexes {ji, U). This is indeed a natural question, quoted as Research 
Problem 2.3.1 in Furthermore, it can be shown, using the stochastic calculus methods in¬ 
troduced in H], as well as the large deviations estimates for semi-martingales given by ||], that 
a moderate deviation principle could be obtained for Alog(Zjv) as long as we can prove the 
boundedness in N of some exponential moments for the random variable 
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This, in turn, can be easily reduced to the asymptotic behavior of quantitites of the form 


' k 

j'=i 


for any fc > 1. This article can be understood then as a preliminary step towards the moderate 
deviations principle for Ziq. 

Let us now present the main results of our paper. The general expansion in terms of N for ( 0 ) 
is given by the following theorem: 


Theorem 1.1 For any n > and f3 < 1, there exist constants C{l3,j) such that 

k 




2=1 


C{P,j) 

N3 


O 


1 


_/Vn+l 




Remark 1.2 One would also expect to get some terms of order N ^ in that expansion. Our 
result shows that the coefficients corresponding to those terms are all vanishing. 

As a particular case of our general result, we have the following 

Theorem 1.3 For any n >1 and P < 1, there exist constants C{P,j) such that 


'{Rh) = E 


i=i 


C{P,j) 

Ni 


O 


1 


_/V"+i 


where 


C{PA) = 


(1-/32)’ 


2 ) ^ -f(b+f;),C(/3, 3 ) = if' + 


(l-/32)4 


(l-/32)7 


In fact, Theorem 1.1 will be an easy consequence of the next proposition, whose proof will be 
detailed in the sequel of the paper. 


Proposition 1.4 For any n > and /3 < 1, a term can be written as a 


finite sum of terms of the type 


i=l 




2=1 


with j,r > 0, j + ^ ^ plus a remainder of order O ^ y i ^, where K{P,I) denotes a constant 

depending on /3 and the family of indexes I = {{ji,li),l < i < r}. Moreover, when 3 + § = f then 
j > 0. 

Remark 1.5 Notice that when r = 0 we get the terms of the form with j G N and j > %. 


These terms are the unique ones that appear in the decomposition given in Theorem 1.1 
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Note that the bulk of the computations one has to perform in order to obtain the expansion 
of y{R\_ 2 ) does not allow us to give a closed formula for all the constants C{P,j). However, the 
analysis of the mechanism needed to obtain those fomulae lead us to write a computer program 
performing this kind of calculations. The listing of the program will be given in an Appendix. 

Our paper will be divided as follows: in the next Section, we will recall some basic facts and 
methods concerning the SK model. At Section 3, for sake of readability, we will explain our 
methods by calculating the first two terms of t he expression of v{R\ 2 )- Then, at Section 4, we 
will prove our main result (i.e. Proposition 1.4). Section 5 will be devoted to the explanation of 
the algorithm of our program, which will be given in the Appendix. 

Troughout the paper all constants appearing in the proofs are called K,K{f3),K{P,k) and 
K{f3, k,n), although they may vary from one part to another. We will omitt the dependence on 
the family of indexes I. 


2 The cavity method 


Let us recall now the basic features of the cavity method (taken from [^), that allows to reduce 
integrals with respect to a 7V-spins system to a iV — 1-spins one: set 


P- 



and (.)_ the averaging with respect to Gn at temperature /3_. Let 
n independent copies of a V-spins configuration. For any j £ {1,.. 
= {fP ,e^), where fP G Sw-i and G { —!> !}• Set 


( 2 ) 

n > 1 and cr^,..., cr" be 
.,n}, we will also denote 


9i = gi,N 

P 


9 ip) = 


Ari/2 


E 


Qi^i’ 


i<N-l 

Then it is easily shown that for one configuration ct, setting crjv = £, we have 

-PHNicr) = -P-HN-iip) + £9 ip ), 

and hence, for / : R, we can decompose (/) into 

(Av/£:„}_ 

\J) ' 

where, for a function <^3 : { —1,1}" ^ R, we set 

1 


Av(^ = 




S’! • • )£n — IT}” 


and 


En = exp 


. l<.n 


Z = (Av£:i)_ = (cosh(g(p')))_ 
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The path between the original spins system and a system in which the last spin is decoupled 
from the other ones, can be defined as follows: for t € [0,1], set 

fft(p) = 

For n > 1, / : ^ R and t £ [0,1], let then (/)* be the quantity defined by 

if), = zr (Av/£:„,t)_ , 

where 


£n,t = exp 


, l<ri 


Zt = {Av£i^t)_ = {cosh {gt{p‘))))_. 

Set also vtif) = E[{f)t\ for t G [0,1], and notice that v{f) = vi{f). 

The following propositions, that we label for further use, are then shown in 

Proposition 2.1 For a function f on and t>0, we have 

v'tU) = 0^ MfsFi'Rpi') 

<n 

,n+l) 

l<n 

2n(n+ 1) 


+ 


~^t{f ^n+l^n+2-^n+l,n+2) ■ 


Proposition 2.2 If f is a function on E^, then 


N2 


3 Expansion of h'{Ri 2 ) 

We begin by introducing some notation. 

Definition 1 Given an expression Y[i=i ^ji^uRji,! define the e-order of this expression as the 
number of different indexes that appears an odd number of times in the set {jifii; 1 < i < k}. 
Notice that the e-order can be computed as 

k 

mod + 

n i—1 

Notice that the e-order is always a nonnegative even number. Note also that 

k k 

Y\_^jNhRji,h and 

i=l i=l 
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are of the same £-order. 

In order to get the expansion of v{R\ 2 ) we need the following technical lemma that allows us 
to obtain of a lot of terms. 


k 

Lemma 3.1 Given set A its corresponding e-order. Then for all P < 1 

i=l 


(a) 


(b) 

(c) 


Sj-Si-Rj. ^J = 0 when A > 0; 


k 


\MY[sji£iiRj,,u)\ < when A = 0, 

i=i ^ ^ 




2=1 

k 




2=1 


K{k,P) 

TVi, 

K{k,P) 

Ni. 


Proof: The case (a) when A > 0 follows immediately from the definition of vq. 

The other cases follow easily from the estimates v{Ri^ 2 ) — ; ^((^ 1 ^ 2 )^^) ^ and 

vt{{Rf^ 2 Y^) — P)'^{{Ri. 2 YY given by Talagrand (see (2.88), (2.89) and (2.91) in [^), and 

the following version of Holder’s inequality: if ti, r 2 > 1,^ ^ = 1, then for all t G [0,1] 

vt{hf2) < ■ 


□ 


Then we are able to prove the following expansion. 
Proposition 3.2 For all P <1 we have 

1 P^{l + P^) 


^(^ 1 , 2 ) — 


O 


1 

Ni 


Nil-P"^) iV2(l-/32)4 
Proof: Using the symmetry property and the relation eie 2 i?i ,2 = ;^ + ^i^ 2 -Rr 2 

^{Ria) = ^{.£i£2Ri, 2) = -^ + ’^{£l£2Rf2)^ 


( 3 ) 


Performing a Taylor expansion and using Proposition ^ and Lemma |3.1| in order to bound the 
last term we get 


h'{eie2Ri 2) — ^'o(eie2l?i^2) +'2o(eie2i?i 2 ) 

+ 2^"o(£ie2l?i,2) + q^o^\^i^2Ri^2) 


1 

TTs" 

N 2 
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( 4 ) 










Notice first that ^'o(ei£ 2 ^i 2 ) ~ 0. To deal with the fc-derivative term we will apply Proposition 


2.1 k times iteratively and we simplify all the null terms using the first part of Lemma 3.1 
Study of u'o{eie 2 Ri 2 )- It is easy to check that 

iy'o{eiS2RZ2) = P‘^'^oiiRf2)^)- 

Using again a Taylor’s expansion we can write 

iy'o{eie2Rf2) = 


= 


’^{{Rf2r) - + O 


iVi 


( 5 ) 


From Proposition 2.1 and Lemma 3.1 we get that i^'o{{R^ 2 )^) = 0- On the other hand, using again 
the relation ei£ 2 l?i ,2 = ;^ + Si^ 2 Ri 2 ^ ‘^nn write 


^{(.^■1.2?) = HR 1 . 2 ) - j;j’^{£i£2Ri,2) + 


Ar 2 - 


( 6 ) 


Finally 


^"om.2f) = h((i?r:2)')-4^^o((i?i:2)"(^r:3)")+3i^o((i?r:2)"(i?3:4)")] 

= /?' [^( 1 ^ 1 , 2 ) - MRhRh) + MRhRh)] + o , ( 7 ) 

where the first equality follows from the computation of the derivative and the second one is again 
an easy consequence of the relation £i£ 277 i ,2 = + eie 2 Rf 2 - Then, putting together (|), d) and 

and using t 2 {eie 2 Ri, 2 ) = v{Ri. 2 )^ we obtain 


0/42 o2 

Meie2Rl2) = PMRI 2 ) - -^’^(Rh) + ^ 


N 


iV 2 


-Ip^HRU) + - Ip^i^iRhRh) + O . ( 8 ) 

Study of v"o{eie 2 Rf 2 )- Using similar arguments, that is, computation of the derivatives, 
Taylor’s expansions and the relation £ie 2 i?i ,2 = ;^ + £i£ 277 i ”2 we can easily prove that 

12 o(£i£27?i^ 2) = —4/3 27?l,3l?2,3) 

= —4/3 [^ 2 ( 7 ?! 27^l,3l?2,3) ~ ^o(I?i^2l^l,3l^2,3)] 4" O 

+ 0 

1 


= -4/3"* v{Ri^2Ri, 3R2,3) — j;^^{Rl,2Rl,3£2£3) 


1 

Ni 


= -WuiRl,2Rl,3R2,3) + O 


Ni 


( 9 ) 


Notice only that we have used that voiRi 2 R 1 3 R 2 3 ) = 0 by Proposition 2.2 and Lemma 3.1 and 
that by the symmetry property and again by Lemma 3.1 


|i^(7?i,27?i,3£2£3)| — |j^(7?i,27?i,37?2,3)| < — T- 

N 2 
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Study of Uq {eiS2Rf2)- similar arguments we can now obtain 


^o\si£2Ri^2) ~ [^o{{Ri^2) ) 12 r ' o (( i ? 2 , 2 ) (^2,3) ) 

+9r'o((i?i 2 ) (^ 3 , 4 ) ) + 36r'o(i?]^ 2^1,3^2,4'^3,4)] 

= HRt2) - 12,yiRl^Rl,) + MR^Rli) 


+36l2(i?i 2i?l,3^2,4-R3,4)] + O ( ) . 


(10) 


Using the same ideas, it is also easy to check that |i^(i?i, 2 -Ri, 3 ^ 2 , 4 ^ 3 , 4 )| < Then, putting 




together (|), (|), (|), (|) and (|I^, and observing that some terms cancels out, we get 

1 0/32 

HRh) = j;f+P^'^iRl2) + ^-^Heie2Ri,2) 


7V2 N 


or the equivalent equality 

* 2 (^ 1 , 2 ) = 


-2(3MRi,2Ru 3R2,3) - l^'^’^iRU) + O 


\N^ N'^ N ^ 


1-/3' 


-2pMR,aRi3R2,3) - ^/3®^^«2)) + O . 

Study of i'{Rf 2 )- Applying the same type of arguments it is easy to get 
^{Ria) = J^(ei£ 2 l??, 2 ) 

= J^(£ie2(l?i,2)^) + jif^iiRia')'^) + C> 

= r'o(er£2(l?i,2)^) + ^o(£re2(l?i,2)^) + jq^{{Ri,2f) 

= /3^t'o(l?i,2) + 2)) + 13 


( 11 ) 


iVi 


— P^^o(Ri,2) + 


+ 0 


TvrS J ■> 

N 2 


A^2(l _/32) 

where we have used that implies that v{R\ 2 ) = A O il is clear that 


K< 2 ) = 


+ 0 


TvrS } ■) 

N 2 


/V 2 (l -/32)2 

Study of v{Ri^ 2 Ri, 3 R 2 , 3 )- Using the same ideas it is not difficult to obtain 


y[Ri.2Ri,3R2,3) = 


A^ 2 (l _/32)3 


+ 0 


3 . r 5 / ? 

N 2 f 


( 12 ) 


( 13 ) 
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Then, plugging (|^ and (|^) into and applying iteratively this new equality we get 

1 /32 2/32 1 




iV(l-/32) 7V2 (i_^2) tv iV(l-/32)2 


2/34 


/3® 


3V2(1_/32)4 3 TV2(1-/32)3 

1 /32 + /34 „ / 1 


O 


1 

TVi 


3V(l-/32) 3V2 (i_t32)4 


O 


ATS 


□ 


4 Proof of the main result 


In this section we give the proof of Proposition 1.4. This proof is based on a sequence of trans¬ 
formations that generalizes the tools used in the proof of the expansion of v{R\ 2 )- We begin by 
introducing a new dehnition and explaining the four transformations we will use. 

Definition 2 We will say that a term is of order m if its absolute value can be bounded by a 
constant times . 

^ fc 1 ^ 

Remark 4.1 From Lemma 3.1, notice that the terms [ Sj^ei^Rf j.) 

/V J .X. .X. /v j X. X. Ji 1 1 


2=1 


2=1 


and —vo(Y\_£ji£hRji,ii) are of order j + |. 


Let us now introduce the four transformations that we will use througout this paper. 


Transformation a. From R to R . 

Using the relation £ji£uRji,ii = SjFuRZ i- F get 


-Ji.ii ' AT’ 
k , k 


(YY^ji^uRjiM) — -ffk + 


TVfc 


where is a finite sum of terms of the type 


k—m 






i=l 


with fc — 1 > TO > 1, that are of order greater than or equal to (see Remark 4.1). 

Transformation b. Taylor’s expansion. 

For a fixed n > using Taylor’s formula, we obtain 


2n — k 


1 




2=1 


2=1 


2=1 
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Applying Proposition ^ m times we get that ^ finite sum of terms of 

the form 

k m 


i' = l 


where I' designates the sequence {{ji,li), 1 < * < m}. 
So we can write 


— ^oiWsji^hRji^u) + ^'o(W_£3i^hRji,u) + ?b + O 


/ 1 


XN'^+h 


i—1 i—1 i—1 

where Tf, is a finite sum of terms of the type 

K{(i)vo(Wej,euR-i^) 

i^l 

with 2n > p > k + 2. These terms are all of order greater than or equal to | + 1 (see Remark 0 ). 

Transformation c. Inverse Taylor’s expansion. 

Fixing again n > |, using a Taylor’s formula we get 


2n—k ^ k 


^o( 


i—1 i—1 m—1 ' i—1 

Using again the expression we can write 

k k 


iV "+2 


2=1 


2=1 
2n—k 




m—1 I' 


2=1 

k 


i=i 


Repeating this procedure with all the terms r/olY^Sj^ei-R^. i. obtain that 


i=l 


i=i 


Vo 


(Y\_^Oi£hRji,ii) — + Tc + 0 


i—1 i—1 

where is a finite sum of terms of the type 


iV”+! 


K{(i)v{\\_£3i^kRj„i,) 

2=1 

with 2n>p>k + l that are of order greater than or equal to (see Remark 0). 
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Transformation d. From R to R 

Using the relation Sj^ei^RJ^ i, = £ji£iiRji,ii — we get now that 

k 1 ^ 

^iY[^n^hRj,,u) = ( - Trf + + Td, 


2=1 


2=1 


where Td is a finite sum of terms of the type 


k—r. 




2=1 


with /c — 1 > m > 1, that are of order greater than or equal to (see Remark 4.1). 

The following remarks explain how can we apply these transformations. 

Remark 4.2 For any fixed n,j, k integers applying transformation c and transformation d we can 
develop the expression 

1 




into a remainder of order n + ^ plus a finite sum of terms of the type 

^ m 


2=1 


with n > I ^ ^J+ 2 ’ ^ T m > j k, I >j. Notice again that all this new terms are of order 
greater than the initial one, except for the term 

1 ^ 

j^iy([[£n£hRnM) 


2 = 1 


that is of the same order. 


Remark 4.3 Fixed n, j, k, using transformation b and Remark we can expand the expression 

1 1^1 


2 = 1 


into a remainder of order n + ^ plus a finite sum of terms of the type 

- m 

2=1 

amb n > I + ^ + I + m > j + k, I > j. Notice also that all this new terms are of order 

greater than the initial one, except the term — ■ ^(Y\_£ji^iiRji,ii) of initial order. 


i=l 
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1 ^ 

Moreover, if we denote by X the e-order of our initial expression i )i 

2—1 

check that KQ{f3) = ff^ when A = 0 and KQ{f3) = 0 when A > 0. 

k 

Proof: [Proposition Set A for the e-order associated to W_£ji£iiRji,ii- We need only to 

consider 3 cases. 

• Case a): A > 4. 

Applying Transformation a we obtain that 

k k 


2=1 


+T;*, 


i=l 


i=l 


where T* is a finite sum of terms of the type 

1 


N3 


jHY[e3i£hRj,,0 


2=1 


with j ^ Notice that all these terms can be studied by Remark 4.3. 

k 

So, it suffices to deal with the term viY^Sj^e^Rj, i,). Since A > 4 it is easily seen that 


2=1 2=1 

and moreover, applying Transformation h we get a finite sum of terms of the type 

r 

KiP)MY[£ji£hRj„i,) 


with r > k 2. Remark 4.2 completes the proof of this case. 
• Case b): A = 2. 

Applying Transformation a we can write 

k k 


'(Y[^neuR,,,0 = Hl[^nei,R-,u) + K* + 


2=1 


2=1 


where T** is a finite sum of terms of the type 


k-l 


N 




2=1 


(15) 
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and S** is a finite sum of terms of the type 




with j > l,r+j = k, |+J> f + 1. 

k 

Let us study first the term i^)- Since A = 2 a trivial verification shows that 

k 

(H /J = 0 there exists jA+i, lk-\-i verifying 






2=1 




2 = 1 


and such that the e-order of (H is zero. So, applying Trans- 

formation b we can write 

k k 


k'Y[£ji£iiRj,,u) - f^'^M{f[£ji£iiRj,,k{£jk+i£h+iRj^+^,i^+J)+Tr* + 0 ^ 


(16) 


where T^** is a finite sum of terms of the type 

r 

KWMY[£3i£iiRji,ii) 


i=l 


with r > k -\- 2. 


Let us study now the term /3^r'o ((nil eii£ii%,z,)(£>+i£ 4 +i-Rj„+i. 4 +i))- Notice first that 

kiiUki ^ji^iiRJ,,u)i£jk+i£ik+iRJ^+^,i^+J) = 0- Then applying Transformations c and d we 
can write 

k 

P^M{Y[£ji£hRj,,0{£3,+i£h+,Rj,+,u+i^) 

k 

= l3‘^’£i(Y\_£ji£iiRji,ii){£jk+i£h+iRjk+i,ik+i)) + Rk + O ( ^ , ( 17 ) 

2=1 ^ ' 


2=1 


where T*^ is a sum finite of terms of the type 

£h£h Rji m) 

2=1 

with n > j + § > 
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Finally, by a symmetry argument it is clear that 

k k 

^((J_ J_ ^((JT ^3i^i)^3k+i^k+i) 

2=1 2=1 




(18) 


2=1 


2=1 


Putting together (p^, ([l^, ( [1^ ) and ( p^ we obtain 

k k 

+ t:* + S:* 

2 = 1 2=1 

or equivalently 


k 

i=l 


1 

1-/32 



ry-i*** 

■ ^ b 




cd 


o 




Applying Remark 4.3 to the terms of T** and S'**, and Remark 4.2 to the terms of T*** it 


is clear that we only obtain a sum of terms of the type 






with j + L>k±i, 

Finally observe that the terms satisfying J come from T** and so they have j = 1. 

• Case c): A = 0. 

By symmetry we have 

k k—1 

i=l i=l 

The expression 11^=1 ^ji^hRji,k is of e-order equal to 2, so we can reduce this case to the 
case b). Notice that since we have now order applying the case b), we will obtain terms 
of order | , but with j > 0. 


□ 


5 Computer Program 

In the proof given in the above section we found the basis of the algorithm that allows us to 

k 

compute the expansion of a term v(W£ji£iiRji,ii)- 

2=1 

The algorithm is the following 
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• 1. Study of e-order. When the e-order is equal to 0 or to 2 we keep away part of the term 

_^ 

1 -/ 

also use a symmetry argument. 

• 2. From Rj^i to i?“;. (transformation a). 

k k 

• 3. From to ^'o((transformation b) 

1—1 i—1 

k k 

• 4. From i/o(JJq) to i^(JJej.eqi?” (transformation c) 

i=l 1=1 

• 5. From R~^ to Rj^i- (transformation d). 

After each one of these steps we have a procedure that simplifies the terms obtained. That is, 
we eliminate the terms with greater order than the grade of our expansion. In order to do it we 
give a refinement of Lemma [3.l| that allows us to simplify more terms. 

k 

Proposition 5.1 Given ^ corresponding e-order. Then 

2=1 


k 

will bo simplified later. When the e-order is equal to 0 we 


(a) < 


K 


i=l 

k 


.- 7- fc I X 5 

iV2 + 4 


(b) W(Y[ej^euRj,,u)\ < 


K 


-m -r k. I X 

iV2 + 4 


Proof: 


The case A = 0 has been studied in Lemma 3.1. Assume now A > 0. 
Let us study first (a). Using a Taylor’s expansion we get 




2 = 1 


2 = 1 




n—1 2=1 


V2/- i=l 


for some t G [0,1]. Notice that all the terms except the last one have strictly positive e-order. So, 
from Lemma |3.1| we easily get 


V2/- i=l 


1=1 
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From Proposition 2.2 and Lemma 3.1 we get 




(H ^ ^ I ( n r 


2 \ 2 


Z =1 


Ni 


\i^l 
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K f K 


Ni \N^ 
K 


, _ X I fc ' 

A/^4 + 2- 


To deal with (b) we will use the relation Sj■el^Rj^^l^ = SjiShRj. i- + ^- Then 


k k / ^ k-j 


i=l j=0 


I i=l 


k-j 

Notice that all the terms Sj^ei.R~ are of e-order greater than or equal to max(0, A — 2j). So, 

i—1 

using the inequality obtained for (a) we have 


k-j 


Ni 




i=l 


< 


K 


< 


K 


m + ^ ■ ivt+^ 


□ 

In the Appendix you can find the main parts of a computer program where we have developped 
this algorithm. This computer program has been written in Maple?. The program is of free use 

k 

and is available from the authors. One must introduce the expression ^'(Y^ej^ei^Rj^j.) using 

i—1 

the notation [[ji, ^i], [j 2 j ^ 2 ], • ■ •, [jfc) ^fe]]- Then to call the computer program and to obtain the 
expansion of such expression until order m, one has to write Zi], [j 2 , h], ■ ■ ■, [jfe, ^fc]]]) w). 


Remark 5.2 For instance we can get 


i^(e 2 e 2 ??i,2^1,3) 
J^(Ri,2Ri,3-Z?2,3^2,4) 
v{£l£2£j,£iRl,2R3,i) 


1 


(1 -/32)37V2 

1 -f /? 

(1 -/32)5iV3 

1 

(1-/32)27V2 


o 

O 


1 

1 


tot([[[l,2],[l,3]]],2); 

tot([[[l,2],[l,3],[2,4],[3,4]]],3); 

tot([[[l,2],[3,4]]],2). 


Appendix 

Main program (tot) Parameters; expres indicates the expression we want to expand and derivord 
is the order of the expansion 
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> tot:=proc(expres,derivord) 

local taula,coef,i,k,m; 
global resfin,resfin2,kdes,hnm; 

taula := expres; coef := [1]; resfin:=0; resfin2:=0; kdes:=derivord; hnm:=l; 
maxima(taula[l]); k:=0; 
for i from 1 to r do; 

• m;=(numboccur(taula[j],i) mod 2); if m=l then k:=k+l; end if; 
end do; 

if kdes < ((l/2)*nops(taula[l]))+k/2 then hnm:=0; resfin2~0; end if; 
while (nops(taula)*hnm)>0 do 

inici (taula,coef); taula:=tauls; coef;=coefs; resfin:=collect((resfin+resu),l/N); 
end do; 

for i from 1 to kdes do resfin2:=resfin2+(simplify(coeff(resfin,N,-i))*N" (-i)); end do; 
print (RESULT, resfin2); 

end proc; 

Procedure inici This procedure calls the four transformations explained in the paper 

Procedure rseparar We study the value of the e-order associated to each term. When it is equal 
to 0 or to 2 we keep away a part of the term that will be simplified later. 

Procedure maxima Compute the maximum of a list 
Procedure rmesamenys Transformation a 

> rmesamenys:=proc(taull,coefl) 

local i,j,ll,k,l,m,kter,long,coefll,coef,taula,taulab,fer,coefb; 
global taulabb,coefbb,resb; 

taula:=taull; coef;=coefl; resb;=0; taulab~[]; coefb:=[]; 
for i from 1 to nops(taula) do; 
long— nops(taula[i]); 
for j from 0 to 2" long-1 do; 

kter—0; conversio(j,long); lh=[]; coefll—coef[i]; 
for k from 1 to long do; 

if jbase[k]=0 then coefll;=coefll*(l/N); 
else kter:=kter-|-(l/2); ll:=[op(ll),taula[i][k]]; end if; 
end do; 

kter;=kter-degree(coefll,N)-|-(comptar(ll)/4); 
if kter<kdes-|-(l/2) and 11<>[] then 

taulab:=[op(taulab),ll]; coefb;=[op(coefb),coefll]; 
end if; 

if kter<kdes-|-(l/2) and 11=[] then resb:=resb-|-coefll; end if; 
end do; 
end do; 
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taulabb := []; coefbb := []; taulabb := [op(taulabb),taulab[l]]; 
coefbb := [op(coefbb),coefb[l]]; 
for i from 2 to nops(taulab) do 

for j from 1 to nops(taulabb) while (taulab[i]<>taulabb[j] or 
degree(coefb[i],N)<>degree(coefbb[j],N)) do 

if j=nops(taulabb) then taulabb := [op(taulabb), taulab[i]]; 

coefbb := [op(coefbb), 0]; 
end if; 
end do; 

fer:= coefbb[j]+coefb[i]; coefbb:=subsop(j=fer,coefbb); 
end do; 
end proc; 

Procedure conversio Convert a number to a binary expression given in a list 
Procedure comptar It eounts the ocurrences of each number in a list 
Procedure simplif It joins all the terms with the same expression 
Procedure vdesenvola Transformation b 
> vdesenvola~proc(taull,coefl) 

local taula,coef,taula2,taula3,vtaula4,vcoef4,coef2,coef3, 
vtaula3,vcoef2,vcoef3,coef5,coeft,ll,i,l,j,ik,k,x,kh,kter; 
global tauldh, coefdh; 

tauldh:=[]; coefdh:=[]; kh:=kdes; taula:=taull; coef:=coefl; taula2~[]; coef2;=[]; 
for i from 1 to nops(taula) do; 

kter:=2*(kh-((nops(taula[i])/2)-degree(coef[i],N))); 
vtaula3:=<< 0 >>; vcoef3;=<<coef[i]>>; 

for ik from 1 to nops(taula[i]) do; vtaula3;=<<vtaula3 [ taula[i][ik]>>; end do; 
vtaula3:=SubMatrix(vtaula3,l--l,2..nops(taula[i]) + l); 
vtaula4:=vtaula3; vcoef4:=vcoef3; 
vanulareps(vtaula4,vcoef4,l); 

tanldh;=[op(tauldh) ,op(taulde)]; coefdh:=[op(coefdh) ,op(coefde)[; 
for j from 1 to kter do; 
vderivar(vtaula3,vcoef3); 
vtaula3:=vtauld; vcoef3;=vcoefd; 
vanulareps (vtauld, vcoefd J); 

tauldh:=[op(tauldh),op(taulde)]; coefdh:=[op(coefdh),op(coefde)]; 
end do; 
end do; 
end proc; 


Procedure vanulareps It e anc els all the terms that are null or that are of order greater that our 
expansion (using Proposition 5.1) 
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Procedure vderivar It gives us the derivative of an expression 
Procedure mmaxima Compute the maximum of a table 
Procedure vdesenvolb Transformation c 
> vdesenvolb:=proc(taull,coefi) 

local taula,taula2,taula3,taula4,tauldh,coefdh,ik,vtaula4, 
vcoef4,coef,coef2,coef4,coef5,ll,i,l,j,k,x,kj,kter,n; 
global tauldj, coefdj; 

tauldj:=[]; coefdj~[]; taula4:=[]; taula2:=[]; coef2:=[]; coef4:=[]; 
kj:=kdes; taula:=taull; coef;=coefi; 
for n from 1 while nops(taula)>0 do; 
tauldh:=[]; coefdh;=[]; 

tauldj;=[op(tauldj),op(taula)]; coefdj;=[op(coefdj),op(coef)]; 

taula3:=taula; taula:=[]; 

for 1 from 1 to nops(taula3) do; 

kter;=2*(kj-((nops(taula3[i])/2)-degree(coef[i],N))); 
vtaula4;=<< 0 >>; vcoef4:=<<-coef[i]>>; 

for ik from 1 to nops(taula3[i]) do; vtaula4:=<<vtaula4|taula3[i][ik]>>; end do; 
vtaula4—SubMatrix(vtaula4,l..l,2..nops(taula3[i]) + l); 
for j from 1 to kter do; 
vderivar (vtaula4,vcoef4); 
vtaula4:=vtauld; vcoef4:=vcoefd; 
vanulareps(vtauld,vcoefd,j); 

tauldh ~ [op (tauldh) ,op (taulde) ]; coefdh:=[op (coefdh) ,op (coefde) [; 
end do; 

taula4:=[]; coef4:=[]; 
end do; 

tanla~tauldh; coef:=coefdh; 
end do; 
end proc; 

Procedure rmenysames Transformation d 
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